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Abstract. The fragmentation of a two-dimensional circular disc by lateral impact is investigated using 
a cell model of brittle solid. The disc is composed of numerous unbreakable randomly shaped convex 
polygons connected together by simple elastic beams that break when bent or stretched beyond a certain 
limit. We found that the fragment mass distribution follows a power law with an exponent close to 2 
independent of the system size. We also observed two types of crack patterns: radial cracks starting from 
the impact point and cracks perpendicular to the radial ones. Simulations revealed that there exists a 
critical projectile energy, above which the target breaks into numerous smaller pieces, and below which 
it suffers only damage in the form of cracks. Our theoretical results are in a reasonable agreement with 
recent experimental findings on the fragmentation of discs. 

PACS. 46.50.+a Fracture mechanics, fatigue and cracks - 62.20.Mk Fatigue, brittleness, fracture, and 
cracks - 64.60.-i General studies of phase transitions 



1 Introduction 



Fragmentation of finite size systems is a widespread phe- 
nomenon in nature which is of considerable scientific and 
industrial interest. Despite the intense research in various 
fields of science and technology, complete analysis and un- 
derstanding of fragmentation has not yet been achieved. 
In industry comminution, or the process of size reduction 
of granular materials, is very important. There is particu- 
lar interest in the net energy required to achieve a certain 
size reduction and the energy distribution of fragments 
during the fragmentation process. For different comminu- 
tion operations like crushing and grinding, early laws were 
developed by Von Rittinger, Kick and Bond 1 , having a 
relatively narrow size range of applicability. Many exper- 
imental and simulational studies have been undertaken 
to extend our understanding of fragmentation to a wider 
range of situations. 

Fragmentation occurs over a wide range, from the col- 
lisional evolution of asteroids [2j to the fragmentation of 
molecules |3I4| . At intermediate scales, there is the degra- 
dation of materials comprising small agglomerates em- 
ployed in process industries such as pharmaceuticals, chem- 
icals, fertilizers and detergents. There are also many geo- 
logical examples associated with the use of explosives for 
mining such as oil shale industry, coal heaps etc. It was 
found in many cases that the probability density function 



F(m) of fragment masses follows a power law 
F(m) ~ m~ T . 



(i) 
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Several numerical 4 5 6I7I8I91 1 01 1 1 1 1 21 1 31 1 41 1 5 j and ex- 
perimental |2I17I18I19I20I21I23I24I25| studies have been 
performed on the fragmentation of solids and particle ag- 
glomerates induced by impact under different conditions. 
In all of these studies one obtains F(m) ~ m~ T for the 
fragment mass distribution with an exponent r depending 
on the dimensionality of the system, i.e. in higher dimen- 
sions usually a larger value of r is obtained j2IHI4l5lh17l8l9ll()ll 1112111111 
Tabled summarizes the most important numerical and ex- 
perimental results on impact fragmentation. 

Recently, "lateral experiments" where a hypervelocity 
projectile strikes the side of a plate was carried out by T. 
Kadono [23]. Moreover, another experiment was carried 
by T. Kadono and M. Arakawa by shooting a cylin- 
drical aluminum projectile with a diameter of 15mm and 
height of 10mm on thin Pyrex glass plate targets. They 
calculated the cumulative number of fragments larger than 
a given mass, and found that the resulting function was 
a power law with an exponent of about 0.60. This func- 
tion is the integral of F(m), so their results correspond to 
r = 1.60. They found two types of crack pattern during 
fragmentation process: radial cracks initiating from the 
impact point and cracks perpendicular to the radial ones. 
They considered that the radial cracks formed prior to the 
perpendicular ones |24l25j . 

In the present paper we have studied the lateral impact 
of a projectile into the right hand side of a disk following 



2 



Bhupalendra Behera et al.: Fragmentation of a circular disc by projectiles 



Table 1. A Summary of recent numerical and experimental 
studies of fragmentation, showing the exponent of the fragment 
mass distribution. 



Numerical Results 


Exponent r 


Explosion of disc-shaped solids 7 10 


2.0 


Impact of a projectile into a block |7I1U| 


1.98 


Collision of two discs 18191 


1.75 - 2.27 


Lennard- Jones solid in d=2 15115 16 


1.4-1.5 


Experimental Results 


Exponent r 


Collisional disruption of ice |17| 


1.64 ±0.06 


Fragmentation of shells |18| 


1.35 


Fragmentation of dry clay plates 




(for large fragments) |19I21| 


1.12 - 1.27 


Fragmentation of dry clay plates 




(for small fragments) |19| 


1.5 - 1.67 


Proton-nucleus collision |2()| 


2.6 


Free fall of glass rods |22I 


1.2-1.8 


Impact on glass tubes fZS\ 


1.5 


Fragmentation of plaster plates |24| 


1.1 - 1.3 


Fragmentation of glass plates 


1.7 


Peripheral Au + Au collisions |28j 


2.2 



the experiments of Refs. 24 25 . Here we have taken a par- 
ticle composed of numerous unbreakable, undeformable, 
randomly shaped polygons which are bonded together by 
elastic beams. This model has been used before for the 
study of the fragmentation process in different physical 
situations |7I8I9I10| . The contacts between the elements 
(polygons) can be broken according to a physical breaking 
rule, which takes into account the stretching and bending 
of the connections. Based on simulations of the model, 
we performed a detailed study of the failure evolution at 
different impact energies and of the nature of the crack 
propagation during the fragmentation process. We also 
study the statistics of average fragment mass and aver- 
age largest mass, the energy and velocity distributions of 
fragments during the fragmentation process when chang- 
ing the imparted energy. 

We have given an outline of the theoretical background 
of the model in Sec. a short description of the basic 
mechanism of fragmentation process and crack propaga- 
tion is presented in Sec. 03 The numerical results at dif- 
ferent initial conditions are given in Sec. 21 where the ve- 
locity, energy and mass distribution of fragments, and the 
transition from damage to fragmentation are studied. 

2 The model 

To study fragmentation of granular solids, we performed 
molecular dynamic simulations (MD) in two dimensions. 
In order to better capture the complex structure of real 
solids, we used arbitrarily shaped convex polygons that 
interact with each other elastically. The model consists of 
three major parts, namely, the construction of a Voronoi 
cellular structure, the introduction of the elastic behavior, 
and finally the breaking of the solid. 

The Voronoi construction, which is a random tessella- 
tion of the plane into convex polygons, was used to reflect 



the grain structure of the solid. One puts a random set of 
points onto the plane and then assigns each point that part 
of the plane which is nearer to it than to any other point. 
In our case, the initial configuration of the polygons was 
constructed using a vectorizable random lattice, which is a 
Voronoi construction with reduced disorder jJHj. One ad- 
vantage of the Voronoi tessellation is that the number of 
neighbors of each polygon is limited which makes the com- 
puter code faster and allows us to simulate larger systems. 
We generate a square block of Voronoi cellular structure 
from which we cut out a disk, special care being taken to 
get a smooth outer surface. With different seed values of 
the Voronoi generator, different samples are obtained with 
differently shaped Voronoi cells. Previously this model has 
been applied to study fragmentation of solids in various 
experimental situations 7 8 9ll()j . 

In the model the polygons are rigid bodies with three 
degrees of freedom in two dimensions: the two coordinates 
of the center of mass and the rotation angle. They are nei- 
ther breakable nor deformable but they can overlap when 
pressed against each other. The overlap represents local 
deformation of the grains. Usually the overlapping poly- 
gons have two intersection points which define the contact 
line. In order to simulate the elastic contact force, we in- 
troduce a repulsive force between touching polygons. This 
force is proportional to the overlapping area A divided 
by a characteristic length L c (-^- = |(^- + ^-), where r\, 

Tj are the radii of circles of the same area of polygons), 
multiplied by a spring constant that is proportional to the 
elastic modulus divided by the characteristic length. The 
direction of the force is chosen to be perpendicular to the 
contact line of the polygons. Further, damping and friction 
of the touching polygons according to Coulomb's friction 
law are also implemented. 

To bond the particles together, the centers of mass of 
neighboring polygons are joined together with beams that 
exert an attractive, restoring force but can break in order 
to model the fragmentation of the solid. The cross section 
of the beams is the length of the common side of the neigh- 
boring polygons and the length is the distance between the 
centers of mass of the two polygons. The Young's modulus 
of the beams and of the particles E p are independent 
of each other. The beams can be broken according to a 
physical breaking rule of the form of the von Mises plas- 
ticity criterion, which takes into account the stretching 
and bending of the connection 

P ,, = (^y +2 HfkW)> 1 , £? >„, (2) 

y ^6, max / "max 

where = At % i jt % i is the longitudinal strain of the beam, 
l and Qi are the rotation angles at the two ends of the 
beams, and €b jrnax and # ma x are the threshold value of the 
two breaking modes. The breaking condition is checked 
at each iteration time step and those beams where condi- 
tion (|2j) holds are removed from the calculation, and never 
restored. The surface of the grains on which beams are 
broken represents cracks. The energy stored in the broken 
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Parameters 


symbol 


Unit 


Value 


Density 


P 


g/cm 6 


5 


Grain Young's modulus 


E p 


dyne/ cm 2 


10 10 


Beams Young's modulus 


E b 


dyne /cm 2 


10 9 


Time step 


dt 


s 


10- 7 


Failure elongation 


£b,max 


% 


3 


Failure bending 


@b, max 


degree 


1 



Table 2. The parameter values used in the simulations. 



beams represents the energy needed to create these new 
crack surfaces inside the solid. 

The time evolution of the system is obtained by nu- 
merically solving Newton's equations of motion of the in- 
dividual polygons (Molecular Dynamics). For the solution 
of the equations we use a Gear Predictor-Corrector scheme 
of fifth order, which means that we have to keep track of 
the coordinates and all their derivatives up to fifth order. 



3 Fragmentation process and Crack 
Propagation 

In the present paper we apply our model to explore the 
properties of the fragmentation process of a circular disc 
due to an impacting projectile. On the right hand edge 
of the disk, one polygon is chosen to be the center of the 
projectile. Then neighboring polygons are detected (shown 
in black in color in the figure). These polygons are all 
given a large initial velocity, directed towards the center 
of the disk, whose magnitude is obtained from the specified 
projectile energy. 

In Fig. E] we present the fragmentation process ob- 
tained by simulation of a disk of radius 20 cm and a pro- 
jectile energy of 1.6 x 10 8 erg. When the disk is struck by 
the projectile, a high compression wave is formed at the 
impact site, where practically all the beams are broken 
and all the fragments are single polygons. As the com- 
pression wave moves towards the center of mass, radial 
cracks form initiating from the impact point (see Fig. [Hi). 
The amplitude of the initial compressive elastic pulse re- 
maining after primary damage strongly depends on the 
initial imparted energy and the amount of primary dam- 
age that occurred in the contact zone. The pulse reflects at 
the free boundary with opposite phase generating an in- 
coming elongation wave (Fig. HJ)). A high stretched zone 
appears near the opposite boundary, which gives rise to 
high breakage of beams and the propagation of cracks ori- 
ented perpendicularly to the radial ones, Fig. [It jJJ. When 
almost all beams are broken in this zone, the incoming 
elongation wave moves further towards the impact point 
while simultaneously extending the cracks perpendicular 
to the radial ones. When the wave reaches the impact 
zone, the remaining energy is converted to kinetic energy 
of the fragments. So, at the end of the simulation many 
single polygons are flying away, back towards the origin 
of the projectile (Fig. Ei). At the end of the simulation 
we collected all the debris and reconstructed the disc to 
investigate the final crack pattern, see Fig. 




(e) 

Fig. 1. Snapshots of the circular solid disc of radius 20 
cm during the fragmentation process at the critical energy 
E c — 1.6 x 10 8 erg. (a) After time 0.00025 sec, the compression 
wave moving inwards creating the radial cracks initiating from 
the impact point, (b) After time 0.00150 sec, the compression 
wave reflected from free boundary surface and both compres- 
sion wave and elongation wave meet each other close to the 
boundary (c) After time 0.00225 sec, a large number of beams 
break, forming cracks perpendicular to the radial ones, (d) Af- 
ter time 0.00540, the shock wave returns to the impact zone, 
(e) Fragments collected at the end of the simulation and put 
back into their initial positions. The intact beams are colored 
according to their longitudinal strain e. 



The final reconstructed state of crack pattern at the 
end of the fragmentation process obtained at different pro- 
jectile energies is shown in Fig. [21 It is observed that at 
very low energy (Fig. |2t) mostly radial cracks are present 
which are the consequence of the expansion of the sample 
perpendicular to the direction of the projectile. In the en- 
ergy range where solely radial cracks occur the disc prac- 
tically keeps its integrity and suffers only damage in the 
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(c) 

Fig. 2. The final reconstructed state of the crack pattern 
at the end of the fragmentation process obtained at different 
projectile energies E Q . (a) E Q = 0.5xl0 8 erg, (b) E Q = 1.3xl0 8 
erg, (c) E = 4.0 x 10 8 erg. 

form of cracks. Transverse cracks appear when the im- 
parted energy approaches a certain critical value E c (Fig. 
Eb). At very high projectile energy (Fig. Et)> a large num- 
ber of radial as well as transverse cracks form and the disc 
breaks into numerous smaller pieces. 

Finally, two kinds of crack patterns are recognized: ra- 
dial cracks initiating from the impact point and cracks 
perpendicular to the radial ones. The radial cracks devel- 
oped prior to the perpendicular ones. The crack propaga- 
tion and the final crack pattern obtained by simulations 
are consistent with the experimental results |24l25j . 



4 Results 

4.1 Transition from damage to fragmentation 

Previous studies on fragmentation identified a transition 
from damage, where the target maintains its identity, to 
fragmentation, where it is broken into many small pieces. 
The existence of these two states in the impact fragmen- 
tation of discs is confirmed by Fig. 03 where the mass of 
the largest and second largest fragment normalized by the 
total mass are plotted as a function of the projectile en- 
ergy E Q . At low energies the largest fragment contains 
nearly all the original mass, and all other fragments are 
much smaller. Thus, there are two classes of fragments: 
the largest mass is the almost intact target, and the other 
fragments are small pieces that were chipped off from the 



target by the projectile in the vicinity of the impact point. 
The target has only been damaged, not destroyed. At high 
energy however, the two largest fragments have compara- 
ble masses and they are much smaller than the total mass 
indicating a complete break up of the disc. In this energy 
range the target has been entirely destroyed by the pro- 
jectile. 

The transition between the two states occurs at a crit- 
ical energy E c , which can be precisely defined as the en- 
ergy that maximizes the second largest fragment. Its value 
was determined E c w 1.6 x 10 8 erg for the system size 
R = 20cm, and E c « 3.3 x 10 8 erg for R = 30cm. It is 
important to note that the curve of the mass of the largest 
fragment has a curvature change from convex to concave, 
the position of which coincides with the maximum of the 
second largest mass, i.e. with the critical energy E c . 

To analyze in more detail the final state of the frag- 
mentation process, we evaluated the so called single event 
moments of fragment masses defined as 

N 

i 

where N denotes the total number of fragments and m max 
is the largest fragment mass. We evaluated the average 
fragment mass M defined as the ratio of the second and 
first moments M = MijM\. In Fig. H|M is plotted as a 
function of the initial projectile energy E Q for system sizes 
R = 20cm and R = 30cm, where each point represents 
an average over 20 simulations. The system again shows 
two distinct regimes, i.e., damaged and fragmented states, 
separated by a sharp maximum of M at the critical energy 
E c . It is important to note that the transition point is 
nearly independent of the system size R when plotted as 
a function of the specific energy, i.e., the ratio of initial 
projectile energy to initial total mass (see Fig. HI a). It can 
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Fig. 4. M 2 /Mi denned by Eq. Q as a function of the spe- 
cific projectile energy (projectile energy divided by the total 
target mass). Two distinct phases can be observed, the dam- 
aged and the fragmented one, separated by a sharp transition. 
The critical energy was obtained as E c = 1.6 x 10 8 erg for the 
system size R = 20cm and E c = 3.3 x 10 8 erg for R = 30cm. 
The transition point is practically independent of system size 
when M2/M1 is plotted as a function of the specific energy. 





10 10"' 
Mass/Total mass 

Fig. 6. The fragment mass histograms for the system size 
R = 30cm with varying initial projectile energy. The straight 
line shows the power law fitted to the curve at the critical 
energy E c = 1.6 X 10 8 erg with an exponent r = 1.93. 



of E c was varied until a straight line was obtained on a 
double logarithmic plot which also provided an indepen- 
dent way of determination of E c from the previous ones. 
Fig. shows that the functional form 



m tot 



\E n — En 



E > E C1 



(4) 



is obtained above the critical point with a value of E c prac- 
tically coinciding with the value determined previously. 
The value of the exponent (3 was obtained as (3 = 0.49 for 
R = 20cm and (3 = 0.45 for R = 30cm. The two values are 
almost the same with a precision of ±0.05 implying that 
(3 is independent of the system size. 



10' 10° 

Initial projectile energy-Critical energy 



10* 



Fig. 5. The largest fragment mass divided by the total mass 
as a function of the distance from the critical point \E Q — E c \ 
for E > E c - The exponent (3 can be obtained as the slope of 
the least square fitted straight line. 



also be observed that increasing the system size makes the 
peak of M sharper, i.e. the width of the peak decreases 
while the height increases, which is typical for continuous 
phase transitions. 

In order to characterize the behavior of the fragment- 
ing system in the vicinity of the critical point we plotted 
the normalized mass of the largest fragment m max /m^ t 
(where m to t is the mass of the target) obtained at dif- 
ferent imparted energies E Q as a function of the distance 
from the critical point \E Q — E c \ for E Q > E c . The value 



4.2 The Fragment Mass Distribution 

The fragment mass histograms F(m) are presented in 
Fig. for the system size R = 30 cm at varying initial 
projectile energy E Q . F is a probability density function 
so that J F(m)dm = 1 its integral is unity. In order to 
resolve the shape of the distributions at a wide range of 
values, logarithmic binning was used, i.e., the binning is 
equidistant on logarithmic scale. The histograms have a 
maximum at small fragment sizes due to the existence 
of single unbreakable polygons. Below the critical energy 
E Q < E c the distributions have a peak at large fragments 
which gradually disappears when approaching the critical 
point due to the break up of large damaged pieces into 
smaller fragments. At the critical point F(m) becomes 
asymptotically a power law over almost three orders of 
magnitude of m. For the two system sizes studied the ex- 
ponents are found to be practically equal and close to 2.0, 
consistent with previous numerical works |7I8I9I1()| . More 
precisely, one obtains r = 1.93 and r = 1.92 for R = 20 
cm and R = 30cm, respectively. Above the critical point 



6 



Bhupalendra Behera et al.: Fragmentation of a circular disc by projectiles 



0.004 



0.003 



> 0.002 



0.001 




-1500 



-1000 -500 500 1000 1500 
Velocity of fragments, Vx (Cm/Sec) 



2000 



0.005 



0.004 



0.003 



0.002 



0.001 




-1500 



-1000 -500 500 1000 

Velocity of fragments, Vy (Cm/Sec) 



1500 



Fig. 7. The distribution of x and y components of the ve- 
locity of the fragments for system size R — 20 cm, for various 
projectile energies. Each curve represents an average over 20 
simulations. 



slight deviations from the power law can be observed in 
the regime of large fragments due to the finite size of the 
system. Since the effect is rather weak no definite form 
of this cut-off function could be deduced, while in other 
types of fragmenting systems an exponential cut-off func- 
tion was obtained |14I15I18) . 



4.3 Velocity distribution 

In Fig.dthe probability distribution functions n(v x ) and 
n(v y ) of the x- and ^/-components of the final fragment 
velocities are plotted. The initial projectile velocity is in 
the negative x direction. One surprising feature of the x- 
component of the velocities is that the majority of frag- 
ments are moving in the positive x direction, i.e., oppo- 
site to the initial direction of the projectile, and there- 
fore opposite to the total momentum of the system. This 
probably occurs because numerous high speed fragments 
are ejected from the impact site, as shown in Fig. [ljl. To 
satisfy momentum conservation, the fewer, more massive 
fragments must be moving in the opposite direction. Thus, 
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Fig. 8. The ratio of the energy dissipated by cracking Ed and 
of the imparted energy E Q as a function of the specific energy 
Eo/mtot. 



mass and velocity are correlated, which is important to re- 
member when we consider the energy distributions. 

The distribution of the ^/-component of the velocity is 
symmetric about v y = as expected from the symmetry 
of the initial conditions. As the energy increases, the dis- 
tribution broadens. It is worth noting that the maximum 
of n(v y ) is largest when the projectile energy is equal to 
E c , which implies that many fragments are formed with 
low kinetic energy at the critical point. Hence, at the crit- 
ical energy the efficiency of fragmentation is highest be- 
cause the largest fraction of the imparted energy is used 
for cracking. Similar behavior is also observed in case of 
n(y x ) except that the maximum moves away from v x = 
when the input energy approaches the critical energy from 
below, and then moves back towards v x = as the pro- 
jectile energy increases because the intensity of the elon- 
gation wave increases gradually beyond the critical point. 



4.4 Fragment energies 

The efficiency of fragmentation can be quantified by the 
ratio of the energy dissipated by cracking Ed to the im- 
parted energy E Q . In Fig. |Hl the ratio Ed/E Q is plotted as 
a function of the specific energy E /m to t- Obviously, the 
number of broken bonds, and hence, Ed is a monotonically 
increasing function of the imparted energy E Q1 however, 
the ratio Ed/ Ed shows a maximum at the critical energy 
E c independent of the system size, see Fig. |H1 The highest 
fraction of E Q used for cracking at the critical point can 
be estimated in the figure to be ^0.42. 

In Fig. we show the energy distribution of the frag- 
ments. Unlike the mass distributions, there is no clear 
power law form. At large energies, it seems that F(E) ~ 
E~ 2 holds, but a closer examination shows that the graphs 
are curved, not straight, in this region. 
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Fig. 9. Energy distribution of fragments at different projectile 
energies for a disc of radius R = 30cm. 

5 Conclusions 

Motivated by recent experimental findings 24 25], we have 
simulated the fragmentation of a disk by lateral impact 
of a projectile. The target was modeled as an assembly 
of unbreakable, undeformable polygons bound together 
by elastic beams. The simulation permits us to follow 
the fragmentation process in detail, allowing us to ex- 
plain the crack patterns observed experimentally 24 25 . 
As the projectile energy increases a transition 8 occurs 
from damage, where the target is not destroyed, to frag- 
mentation, where the target is broken into many pieces. 
The transition clearly shows up in the behavior of the 
largest two fragments and the average fragment mass as 
a function of the imparted energy. We also found that 
the mass distribution F(m) obeys a power law over sev- 
eral orders of magnitude, with an exponent close to 2, in 
agreement with previous numerical results |7I8I9I1()| and 
experimental findings |17I1 9 20 24 25 28] on the fragmen- 
tation of two-dimensional objects. On the other hand, the 
energy distribution shows no such power law. 

As far as the design of comminution processes is con- 
cerned, this work points out the importance of identifying 
the critical energy E c that separates damage and fragmen- 
tation. To maximize efficiency, such processes should be 
designed to deliver energies just above E c . If the energy is 
too low, the target must be struck several times in order 
to achieve a significant size reduction. If the energy is too 
high, much of it will be wasted, as size reduction does not 
increase much above the critical energy. Size reduction is 
usually not done with high speed projectiles, but rather 
through collisions with hard balls or with walls. Therefore, 
consideration of these circumstances is important. 
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